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We study the fate of superconductivity in the vicinity of a class of metallic quantum 
critical points obtained by coupling a Fermi surface to a critical boson. In such 
systems there is a competition between the enhanced pairing tendency due to the 
presence of long-range attractive interactions near criticality, and the suppression 
of superconductivity due to the destruction of the Landau quasiparticles. We show 
that there are regimes in which these two effects offset one another, resulting in a 
novel non-Fermi liquid fixed point with finite, scale invariant, BCS coupling. While 
these interactions lead to substantial superconducting fluctuations, they do not drive 
the system into a superconducting ground state. The metallic quantum critical fixed 
points are connected to the superconducting regime by a continuous phase transition. 
These results are established using a controlled expansion in the deviation from d = 3 
spatial dimensions, as well as in a large number N of internal flavors. We discuss 
the possible relevance of our findings to the phenomenon of superconducting domes 
condensing out of a non-Fermi liquid normal state near quantum critical points. 
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I. INTRODUCTION 

One of the central unresolved issues of modern condensed matter physics involves the 
enhancement of superconductivity near metallic quantum critical point^ 1 5 . Indeed, many 
of the strongly correlated electron materials, such as the cuprate^®, iron pnictide&P, organic^, 
and heavy fermion systems 9 appear to have enhanced superconducting “domes” when they 
are tuned experimentally towards a quantum critical point. And often, the normal state 
exhibits scaling behavior that is inconsistent with Fermi liquid theory. The manner in which 


2 


such “non-Fermi liquid” behavior is related to the enhanced pairing scale remains a long¬ 
standing, yet exciting and actively pursued topic of investigation. 

The reason for the enhancement of pairing near quantum critical points has been known 
for some time: there are induced attractive long-range interactions between electrons near 
the Fermi surface mediated by critical order parameter fluctuations. These interactions are 
long-ranged because of the diverging correlation length at criticality, and, like phonons, or¬ 
der parameter fields mediate attractive forces. However, precisely the same order parameter 
fields have an opposing effect: they tend to destroy the quasiparticles, enhancing their scat¬ 
tering rate relative to their energy. If this second effect is dominant, the effective description 
of such fermion modes is then no longer governed by Fermi liquid theory. The fermion fields 
develop an anomalous dimension, and there is no longer a quasiparticle description of the 
low energy dynamics. The system will then be governed by a non-Fermi liquid fixed point. 
The destruction of the Landau quasiparticle therefore has a pair-breaking effect, weakening 
the superconducting tendency of the system. 

A challenge remains to predict the circumstances under which the enhanced supercon¬ 
ducting interaction dominates, and those in which the fermion anomalous dimension dom¬ 
inates. Furthermore, given the fascinating properties of quantum criticality, it would be 
extremely interesting to construct models of metallic systems exhibiting quantum critical 
points with non-Fermi liquid behavior in the deep IR. However, in most examples so far, the 
superconducting instability sets in before non-Fermi liquid effects become important. As a 
result, the fixed point is fully covered by a superconducting dome, and quantum criticality 
is not observed. 

In this paper we address these questions in a class of quantum metals where the order 
parameter fields condense at zero momentum. Examples include the Ising nematic transi¬ 
tion in metals, which have been argued to be relevant to the phenomenology of both the 
cuprat^ 10 and iron pnictide 11 superconductors. We consider a solvable large N limit where 
exact statements about pairing instabilities can be made. Our first key result is that the 
competition between the long-ranged attraction and the destruction of Landau quasiparti¬ 
cles can lead to a fixed point where the BCS interaction among fermions is finite. This is 
in sharp contrast to the behavior of Fermi liquids, where the BCS coupling flows to zero for 
repulsive forces, or grows indefinitely, leading to a BCS instability if the couplings are attrac¬ 
tive. The finite BCS interaction fixed point here corresponds to a metallic phase with scale 
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invariant interactions in the BCS channel which do not result in Cooper pair condensation. 
As a result, we will be able to exhibit a “naked” fixed point with critical BCS coupling, not 
covered by a superconducting dome. 12 The thermodynamic and transport signatures at such 
finite BCS interaction fixed points are interesting in their own right and may be relevant 
to experiments involving quantum critical metals. We will study these phenomenological 
properties in future work. 

Our second main result is that there is a continuous transition between the regime where 
the RG flows are always towards enhanced superconductivity, and the regime where IR stable 
fixed points with finite BCS couplings occur. Approaching the transition from the critical 
regime, the IR fixed point annihilates against an unstable UV fixed point and disappears. 
From the other side, the superconducting parameter and all its derivatives vanish as we tune 
towards the transition. This is reminiscent of the Berezinski Kosterlitz Thouless (BKT) 
transition, as we will explain in detail below. 

The paper is organized as follows. In section [TT[ we define the bare tree-level action and, 
neglecting superconductivity for the moment, we describe the dominant quantum corrections 
by solving the Schwinger-Dyson (SD) equations for the system. These equations are exact 
to all orders in perturbation theory, for the large N limit we consider. Next, in section 


III we study the system using the renormalization group. We construct a scaling theory 


consistent with the analysis of the SD equations, determine the one loop beta functions, and 


characterize the non-Fermi liquid fixed point. Section IV is devoted to the analysis of the 
BCS interaction: we show that there is a regime where the 4-Fermi coupling flows to a stable 
fixed point, and another regime where it leads to a superconducting instability. We establish 
that both states are connected by a continuous phase transition. Section [V] contains our 
conclusions regarding the phase diagram of the theory and discusses future directions. Some 
technical calculations are presented in two appendices. 


II. EFFECTIVE ACTION AND QUANTUM CORRECTIONS 

In this section we present the classical theory and compute quantum effects. We use the 
Schwinger-Dyson equations for the boson and fermion two point functions and work in a 
large N limit, which will allow us to obtain results that are exact to all orders in perturbation 
theory. Here we focus on the correlation functions that are local on the Fermi surface, while 
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in sections nm and m we take into account the 4-Fermi interaction in the BCS channel. 


A. The model 

In this work we will analyze the quantum theory for a Fermi surface coupled to a gapless 
boson 0. Our starting bare euclidean Lagrangian is 

L = -Tr [(<9 t 0o)" + (V0 O ) 2 ] + 0j o (&t + £o(*V) — fj,p) V’o + + p bcs 

L 0,0 = ^=0o + <?)0oO) 

Lscs = ~W ZI ^ 0i<yP + - M(-p) • ( IL1 ) 

The subscript ‘0’ denotes bare quantities (we will consider the effects of renormalization 
after integrating out high energy modes below). The bare band dispersion of the fermions is 
denoted by £q( k) and for simplicity we consider a rotationally invariant Fermi surface; the 
chemical potential is [If = £o(kp). The sign in Lpcs is such that Ao > 0 corresponds to an 
attractive interaction. 

We will consider a soluble limit of the theory above. First, we introduce an internal 
SU(N) global flavor symmetry (a generalization of spin rotation symmetry) under which 
the fermion fields transform in the fundamental (vector) representation, whereas the bosons 
transform in the adjoint (matrix) representation. We work in the limit N 1 with go and 
Ao fixed; many diagrams will be shown to be subleading, and it will be possible to resum 
exactly the leading quantum corrections. We note that the large N theory here explores a 
distinct asymptotic regime than the standard large N approach to this problem, in which the 
boson remains a scalar while the fermions are fundamental fields of a global flavor symmetry 
group. Furthermore, we work in d = 3 — e spatial dimensions with e 1, namely near the 
critical dimension for the Yukawa coupling. As discussed below, the small parameter e will 
be used to avoid infrared divergences from corrections to the cubic vertex. 

Besides providing limits where quantum corrections simplify, N and e will also affect 
the infrared dynamics of the theory. Our task will be to determine the low energy phase 
diagram of the theory as a function of N and e. Before proceeding to the discussion of 
quantum corrections, let us develop some intuition by comparing the scales of non-Fermi 
liquid effects and superconductivity. Anomalous dimension corrections become important 
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at a scale 


fJ'NFL ~ G 2 7(A) \ 


(II.2) 


where 


7 = 


9 _ « 

24n 2 v ~ 2 


(II.3) 


and g and v are physical couplings (see below). We will deduce this result shortly, but for 
now we just want to explore some of its consequences. On the other hand, the scale of the 
superconducting gap taking into account the enhancement from boson exchange is (see 
and below) 

_ 7T / N 

g 2 V c(A) \ _ 


hs 


(11.4) 


At the quantum critical point described below, a ~ e. We see here the interplay between 
NFL and gap effects: for e <C 1 /N, superconductivity dominates, and the Fermi surface is 
gapped before the NFL regime is reached. However, in the opposite limit N 1/e, /j-nfl 
fi sc and hence we expect (and will find) strong NFL corrections to the superconducting gap. 
This is the range where a new quantum critical point for the BCS interaction will obtain. 
We will show that both regimes are connected by a continuous phase transition that occurs 
when Ne r<«/ 2.. 

Our strategy will be to first determine the dynamics in the non-Fermi liquid regime, cor¬ 
responding to N 1/e. Here we will neglect the superconducting gap, and then check that 
this is a self-consistent approximation. We will then incorporate effects from superconduc¬ 
tivity and will characterize this phase that occurs when 1 C iV « 1/e. This will be done 
by analyzing the RG /3 function for the 4-Fermi BCS coupling. 


B. Quantum corrections 

Let us then begin the analysis of quantum effects focusing on the self-energies and boson- 
fermion coupling. At sufficiently high energies, the scaling behavior of the theory can easily 
be understood. Considering the Yukawa coupling to be small, scaling is constructed about 
the limit wherein a Fermi liquid is decoupled from the order parameter field. Scaling behav¬ 
iors of the fermions and bosons are governed by Fermi liquid theory and Landau-Ginzburg- 
Wilson effective theories respectively. The dimension of any composite operator, such as 
ifi'i/cf) is then immediately known from the decoupled scaling dimensions of these fields. The 
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conclusion of such an analysis is that the Yukawa coupling has a bare scaling dimension 

3 — d e 
~ ~ 2 ' 


fo] = 


(II.5) 


Our next task is to determine the way in which quantum corrections alter this behavior. 

As discussed in many previous works 14 there are many subtleties involved with taking 
a strict large N limit in this class of theories. Indeed, the scaling behavior in the large N 
limit depends quite strongly on the order in which the N —» oo and oo —> 0 limits are taken. 
If the N — y oo limit is taken first, the resulting fixed points obtained only govern behavior at 
intermediate energy scales^, because they only take into account a subset of the important 
quantum corrections. In particular, in this limit, key 0{1/N) quantum corrections that 
qualitatively alter the IR behavior are neglected; these quantum corrections act as effective 
relevant coupling constants that destabilize potential N = oo fixed points. Here, we wish 
to avoid such peculiarities, and build in all the important quantum corrections, even those 
that formally are 1/N corrections, into our theory. 

Our strategy for obtaining a scaling theory will be to look for effects that are exact to 
all orders in perturbation theory at large N. This will be obtained by investigating the 
Schwinger-Dyson (SD) equations for this system. A key simplification of large N is that 
quantum corrections to the cubic vertex are suppressed by an extra power of 1/N compared 
to the tree level term, at fixed go. This allows us to neglect vertex corrections.- 28 In this case 
we find a closed system of SD equations for the boson and fermion self-energies, expressed 
as follows (see Fig. [Tj): 


9o 


H {q 0 ,q)=D (q 0 ,q) - D 0 (q 0 ,q) = — 

E(po,p) = G ,_ 1 (po,p) - G^ipfrp) = -gl 


dkod d k 
(I 2n) d+1 
f dkod d k 


G(ko, k)G(k 0 + qo,k + q) 
G(k 0: k)D(p 0 - k 0l p - k) 


(IL6) 


j (27r) d+1 

where D{G) refer to the exact boson (fermion) propagator. 

There are many possible solutions to these equations, depending on the parameter that 
is being held fixed while taking the large N limit. We will solve these equations, in d = 3 — e 
spatial dimensions, and in the limit where both N — » oo, kp —y oo holding fixed the following 
quantity 

l.d-1 2 

< n - 7 ) 

where q is a constant that depends on the dimensionality of space. The factor k d jG l /vq 
is proportional to the density of states at the Fermi energy. Physically, this is the scale 
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wv wv 


FIG. 1. Schwinger-Dyson equations for the fermion (full line) and boson (wiggly line) two-point 
functions in the large N theory. 


below which Landau damping of the boson becomes very important and is equivalent to 
the “Debye” screening scale, below which long range Coulomb interactions are screened in a 
metal. As we will see below, by holding M|, fixed, the theory simplifies substantially in the 
IR. These equations build in the dominant quantum corrections in the large N limit above 
the scale of the superconducting gap. 


Above the scale of M^, there are logarithmic corrections to the fermion self-energy, which 
are both frequency and momentum-dependent. They produce a small anomalous dimension 
(proportional to e) and cause a slight reduction of the Fermi velocity. Both of these effects 
can be seen directly in perturbation theory in the UV. For large enough Mjj, however, these 
effects are subdominant in relation to the Landau damping of the bosonic order parameter 
fields, which is a UV finite, non-local quantum correction to D(q), and therefore is invisible 
in a Wilsonian treatment of the problem. Nevertheless, it strongly affects the dynamic 
scaling of the boson. This quantum correction comes directly from the first equation in 


(II.6) and can be interpreted as coming from resumming the geometric series of fermion 
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bubbles (analogous to RPA): 


D(q 0 ,q) 1 = q 2 0 +q 2 + U(q 0 ,q) . 

The behavior in the regime |g 0 | Vo\q |, which will be relevant for us, is 24 

n(«,,?)» Ml M . 

v\q\ 

The boson propagator then takes the approximate form 

D(q o,?) -1 ~q‘ + Ml 


2 , \ f2 l*?o| 


I 


(II.8) 


(II.9) 


( 11 . 10 ) 


and is characterized by a — 3 scaling q 0 ~ \q | 3 /M|,. M D therefore acts as a crossover 
scale, separating the z = 1 UV behavior from the z = 3 IR behavior in this regime. Our 
analysis will focus on scales smaller than M D . 

Below the scale Mn , the self-energy ceases to have substantial momentum dependence 
and depends mainly on frequency. Therefore, at energies much less than Mp the velocity 
renormalization is primarily due to the fermion anomalous dimension (i.e. the velocity will 
vanish at the fixed point where Z vanishes). The non-zero anomalous dimension will also 
cause the renormalization of the Yukawa coupling (note again that vertex corrections are 
suppressed in our large N limit; only field rescaling due to the anomalous dimension causes 


renormalization of the bare Yukawa coupling). Given (11.10), the solution to the fermion 


SD equation in (II.6) is 


£(po) = -iPo — { m d\Po\) 


-e/3 


„ _ % 
« 0 = 


( 11 . 11 ) 


e ' -- - 127r 2 r> 0 

Recalling that the bare coupling «o has engineering dimension e, the self-energy has en¬ 
gineering dimension 1 (as it should), but the scaling dimension with p 0 is 1 — e/3. We 
emphasize again that this scaling applies below Mp, which plays the role of a UV scale in 
our effective theory. 

We conclude that the effective (quantum) Lagrangian involving the fermion kinetic energy 
and the Yukawa coupling is 


Lefw = —0o [iZ(p 0 )po ~ v 0 p±] 0o + ^=05(A: + g)0o(&)0o(<?) 


9o f. 


where Z(p 0 ) is the quasiparticle residue, 

£(po) 


Z(p 0 ) = 1 




Vn 


= 1 + ^5 (Ml\p„\) 


-</3 


( 11 . 12 ) 


(11.13) 
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The fermion momentum here is decomposed radially towards the Fermi surface, 


p = n{k F +p±), 


(11.14) 


and n is a unit vector that defines the position on the Fermi surface. 


The second term in (11.13) represents the effects of quantum corrections. Below a scale 
Pnfl defined by 


Pnfl — 


3CKr 


3/e 


M 


D i 


(11.15) 


the quasiparticle residue Z is dominated by quantum corrections, and the frequency depen¬ 
dence of the fermion kinetic term has the following behavior: 


Po < Pnfl : Z(p 0 )p 0 » —M D 2e/3 pl e/3 


(11.16) 


This is of the form p\ 2l : in other words, at low energies the fermions develop an anomalous 
dimension, 


2y - 
3 


and the quasiparticle residue becomes 




«/3 


(11.17) 


(11.18) 


C. Physical quantities and IR dynamics 


The bare quantities written above are not physically measurable. Once quantum correc¬ 
tions arise, the physically observable fields and couplings depend on the energy scale. In 
the large N limit this occurs only because of 1) Landau damping of the bosons and 2) a 
non-zero fermion anomalous dimension. 

The physical quantities are obtained after canonically normalizing the fermion fields at a 
scale p,, if) = 0 . The physically observable Fermi velocity (as seen in heat capacity 

or in tunneling density of states, for instance) will depend on the energy scale /i as 


v(p) = Z (fi)v Q 


(11.19) 


Note that for /i <C Pnfl, the velocity flows to zero with a power-law determined by (11.18), 

e/3 


v(p) ~ v 0 

\PnflJ 


( 11 . 20 ) 
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Next, we have to determine the relation between the physical dimensionless coupling g(g) 
and g 0 or, more uscfullyp” 

n‘ 2 (n\ 

( 11 . 21 ) 


/ \ _ 9 2 {g) 

a(g) = 


127 f 2 v(g) 

From the effective Lagrangian in terms of canonical fields, we read off a(g) oc Z(/i)~ 1 a 0 . 


Here two powers of Z(g) 1 come from g 2 (g) in (11.21), while v(g) gives an additional factor 


of Z(g) according to (11.19). 

It remains to determine the dimension 5 of a(g). For this, we write the relation between 
bare and renormalized couplings as 


a 0 = M e D 5 Z(g) g b a(g ), 


( 11 . 22 ) 


where Mo is added to match engineering dimensions. The dimension 5 by definition cancels 
factors of the external frequency in loop integrals, such that we get a perturbative expansion 
in terms of a and dimensionless ratios g/po- The one loop fermion self-energy calculated 
using renormalized perturbation theory is 


yv \ • 3 cx(g) --e/ 3-5 g 

£(po) = -Wo - Mr 


L D 


(11.23) 


e " bo| e/3 ’ 

Therefore, 5 = e/3, namely 

Hg)} = | (H.24) 

and, as expected, the microscopic scale M D cancels when working in terms of the physical 
coupling. 


Combining this with (11.22) obtains the relation between «o and a(g) to all-orders in the 
large N theory, 


a(g) = 


«o 


(11.25) 


+ ( M ^) £/3 ' 

This energy dependence of a(g) has the property that for scales g gNFL the coupling 
flows to 

g > (11.26) 


a non-Fermi liquid fixed point that will be analyzed in more detail in Sec. Ill 


Eqs. (11.19) and (11.25) are the main results of this section. They tell us how the physical 


fermion velocity and Yukawa coupling at a scale g behaves. Thus, while perturbation theory 
about the decoupled fermion-boson limit would suggest that [a] = e, quantum corrections 
treated by the SD equations instruct us how these classical dimensions are altered. 
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To summarize: using the SD equations we have computed the quantum corrections that 
occur to all orders in perturbation theory in the large N limit, at scales E < Mp, and 
with M d fixed. Our next objective is to obtain this behavior using scaling and Wilsonian 
RG below the Debye scale. This method will then allow us to determine the evolution of 
the BCS interaction. An alternative analysis of the superconducting gap in terms of its 
Schwinger-Dyson equation (the Eliashberg equation) will be presented in [30], 

III. RENORMALIZATION GROUP APPROACH INCLUDING BCS 

COUPLINGS 

We next determine the RG j3 functions that define the flow of physical couplings for our 
theory. The first step will be to reproduce the SD results in a Wilsonian RG framework. 
We will then focus on the RG flow for the BCS coupling. A Schwinger-Dyson treatment of 
superconductivity effects including the anomalous dimension is more involved, and will be 
presented in 30 . Before proceeding, we note that in the theory in d = 3 — e spatial dimensions, 
it will be computationally more convenient to organize quantum corrections in powers of 
1/e (effectively using e as a regulator), instead of employing the physical cutoffs A/ and A&. 
The map between both approaches was given in 25 . 


A. Scaling theory 


The first step in an RG approach is to construct a consistent scaling. The solution of 
the SD equations places strong constraints on a scaling theory which, in particular, has 


to reproduce (11.24). We now present the scaling that agrees precisely with the form of 


quantum corrections obtained above. 

Ref. [25] showed that a consistent renormalization of the Fermi surface coupled to a 
massless boson requires two independent decimation procedures: the Fermi surface high 
momentum modes are integrated on shells 


Af — dAf < \p±\ < Af , (III. 1) 

where the fermion momentum is decomposed radially towards the Fermi surface, 

P = fi(k F +P±), (III.2) 
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and h is a unit vector that defines the position on the Fermi surface. On the other hand, 
the boson momentum q is decimated towards the origin, with an independent cutoff A b : 

A b -dA b < \q\ < A b . (III-3) 


These two independent momentum-shell integrations are needed to capture the leading quan¬ 
tum corrections to correlation functions. The reason is that some contributions that look IR 
from the point of view of the fermions, actually come from UV bosonic modes, and hence 
have to be taken into account in the Wilsonian RG. Important consequences of this were 
the running Fermi velocity, and tree-level logarithmic running of 4-Fermi couplings. 

This approach has to be modified if the boson has a nontrivial dynamical exponent z b . 
This is discussed in detail in Appendix [A] First, we find that the scaling of the fermions is 
not changed by z b , 

M = H = l, Mp)] = -3/2. (HI-4) 

The dimension of the 4-Fermi BCS coupling is classically marginal in any dimension, since 
we work in the spherical RG for the fermions. 

On the other hand, the bosonic scaling is modified as follows. Given a patch with an¬ 
gular position n on the Fermi surface, we decompose the boson momentum into orthogonal 
components 

q = hq ± + q\\ . (III. 5) 

We show in the Appendix that the correct scaling obeys 

[go] = [q±] = 1 , [q\\] = 1/3 , [0(g)] = ( m -6) 

for d = 3 — e. 

With these scalings, the classical dimension of g becomes 


W = «- 


(HI. 7) 


Therefore, the dimension of g is nearly marginal with this scaling for the overdamped boson, 
as is also the case above the Landau damping scale. This scaling reproduces the quantum 


result (11.24) obtained from solving the SD equations. 

The near-marginality of g after including Landau damping is important, as it is consis¬ 
tent with a smooth crossover between the undamped and overdamped regimes (see alscf 24 -]), 
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and perturbation theory does not break down. Had we scaled the boson momentum homo¬ 
geneously, the result would have been [g\ = (2 + e )/6 ~ 1/3, giving an order one relevant 
interaction in the overdamped regime. Such a relevant coupling would be inconsistent with 
the results of the SD equations. 


B. Local non-Fermi liquid behavior 

Let us discuss first the RG for correlations that are local on the Fermi surface -the self¬ 
energies and the boson-fermion coupling. The one loop corrections are shown in Fig. [2j 



FIG. 2. One loop corrections to the boson self-energy, fermion self-energy and cubic Yukawa vertex. 

We argued in Sec. [TT] that below the Landau damping scale the boson flows to a Zf, = 3 
dynamical exponent. This is not a Wilsonian effect: it comes from integrating particle-hole 
pairs at the Fermi surface; it is a finite renormalization effect that becomes relevant below 
the scale Me- Here we start from this dressed boson propagator and focus on E < Me, but 
we note that the crossover between zj, = 1 and zj, = 3 is smooth 24 . 

In the Zb = 3 regime, the fermion self-energy depends only on frequency. This results in 
a velocity flowing to zero with a rate determined by the anomalous dimension, Z v = Z~ x . 
Furthermore, quantum corrections to the cubic vertex are suppressed by 1/N at large N, 
and can be neglected in our perturbative framework. Therefore, the one loop /3 functions 
characterizing quadratic and cubic correlators are determined purely in terms of the fermion 
anomalous dimension 7 . Using renormalized perturbation theory at one loop, the anomalous 
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dimension is 21 


7 = 


9 


24:ir 2 v 


(III.8) 


In terms of a defined in (11.21), the one loop f3 functions on a local patch of the Fermi 
surface then become 


dlogZ 

27 = — fi — -= a 

d/i 

d log v 

fi —— = 27 
dfi 

da e 

fi— = --a + 270 . 
d/i S 


(III.9) 


This agrees with the (3 functions obtained from the SD analysis by requiring that the bare 


parameter in the relation (11.25) be independent of the RG scale /i. Note that the one loop 
approximation here is exact at large N. 

This system admits a non-Fermi liquid fixed point, 


a* = x , 7* = x 
6 0 


(III.10) 


that is perturbatively controlled at small e and large N. Although we can make no definite 
prediction for e ~ 1, we note that as e —>■ 1 (i.e. for d = 2 spatial dimensions) this fixed 
point approaches the strongly coupled non-Fermi liquid of 17 . It would be interesting to use 
the e expansion to understand more systematically the e —> 1 limit. 


C. BCS f3 function 

We now want to include the effects from the BCS 4-Fermi coupling, which is classically 
marginal and can destabilize this fixed point. The renormalization of the BCS interaction 
proceeds in two steps 25 31 . First, at tree level the boson exchange gives rise to running of 
the BCS couplings in the angular momentum basis; see Fig. [3j At one loop there are two 
additional contributions: from the anomalous dimension and the BCS one loop diagram. 
This is shown in Fig. |4} 

The one loop beta function then becomes (see alscf 13 -) 

d\ A 2 

% = + < IIL11 > 
(We recall that A > 0 corresponds to an attractive interaction.) To our knowledge, the 
term linear in A has not been included in previous works. It captures the non-Fermi liquid 
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FIG. 3. Tree-level running of the BCS interaction. 





FIG. 4. Diagrammatic one loop contributions to the BCS coupling. 


corrections to fermion scattering on antipodal points of the Fermi surface and hence the 
formation and condensation of Cooper pairs. We will find that it has dramatic effects on 


the IR phase structure of the theory, which we consider in Sec. IV 


There are two and higher loop quantum corrections to the 4-Fermi coupling. At the same 


order in N as in (III. 11) there is a geometric series of fermion bubbles, as well as anomalous 


dimension insertions in internal fermion lines. These are automatically resummed into the 
solution of the RG (3 function. Subleading in N effects come from vertex corrections (as 
before) as well as higher loop contributions containing BCS interactions that are not sums 
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of fermion bubbles. This large N suppression is a consequence of the non-planarity of the 


BCS coupling in our theory. We then conclude that (III. 11) is exact at large N. The last 


aspect to understand is whether the 4-Fermi interaction corrects the local non-Fermi liquid 
behavior of Sec. Ill B — All such contributions are again subleading at large N, again due to 
the non-planarity of the BCS interaction (recall that the non-Fermi liquid behavior above 
arises at the planar level). 

Another way of organizing these quantum corrections is to study the Schwinger-Dyson 


equation for the superconducting gap together with (II. 6 ). This approach will be presented 
uP. with conclusions that are consistent with the present renormalization group results. 


IV. QUANTUM CRITICALITY AND FATE OF SUPERCONDUCTIVITY 


In this section, we consider the consequences of Eq. (III. 11). We start by describing, at a 


heuristic level, the various possible fates of the BCS coupling encoded in this equation. We 
then perform a more detailed RG study and discuss the phase diagram of the theory. 


A. Qualitative analysis of the BCS interaction 

For the present analysis, we treat 7 as an independent parameter to exhibit more clearly 
the effects from anomalous dimension corrections (in our theory, 7 = a/2). Furthermore, at 
this heuristic level we will ignore effects from the running of a; these will be incorporated 
below in a more detailed RG treatment. 

First, for simplicity consider the case of a = 7 = 0. In this case, the scalar is effectively 
absent from the theory and we recover the marginally relevant flows of the BCS coupling 
in a Fermi liquid. Attractive interactions grow under the RG whereas repulsive interactions 
weaken. The only difference here, is that BCS couplings are N suppressed because of the 
large N limit we have taken. Thus, there is a limiting case where ordinary Fermi liquid 
theory is captured. 

Consider next the regime where the term proportional to 7 A is subdominant. In this 
case, the effect of the anomalous dimension can be neglected but the term oc a alters the 
qualitative nature of the superconducting instability. There is an exponentially enhanced 
pairing instability, and the pairing scale far exceeds the scale at which the Landau quasi- 
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particles would have been destroyed. The inverse correlation length at which A 1 in the 
IR is of order 


s 


exp 


N IT 
a 2 


A 


(IV. 1) 


This is what one finds in color superconductivity, and is also what has been recently been 
reported in Ref . 33 in the context of 2d quantum criticality. Next, consider the regime in 
which the anomalous dimension plays the most important role. In this case, one can neglect 
the higher order terms of the /3 function, and one finds a stable non-Fermi liquid metal with 
zero BCS couplings. 

Finally, consider the full expression, where all three terms play an important role. Now, 
there are competing, and offsetting effects between the anomalous dimension, and the en¬ 
hanced pairing tendency. Since the /3 function is quadratic in A, depending on its discrimi¬ 
nant there are three different possibilities, illustrated in Fig. [5] 



FIG. 5. BCS /3 function for different values of the anomalous dimension. The two real roots at 
7 > 7c give stable and unstable fixed points. These merge at 7 = r j c and then annihilate; this gives 
rise to a superconducting instability for 7 < 7 c . 

Neglecting the running of a (to be incorporated shortly), we can write the zeroes as 

A± = Ac (l ± Vl-7 c 2 /7 2 ) (IV.2) 

where Ao = 2n 2 N r y, and q c = yja/N. The quantity q c plays the role of a “critical anomalous 
dimension” which separates two regimes. The fixed points are only physical when 7 > q c , 
that is, when the zeroes of the f3 function occur on the real axis. I 11 this case, there is a UV 
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fixed point A+, and an IR fixed point A_. We should emphasize that the critical behavior 
associated to A_ is qualitatively different from a non-Fermi liquid fixed point that is local 
on the Fermi surface, since it affects correlation functions with support on antipodal points 
of the surface. We will discuss further properties of this fixed point in the next section. 

As 7 —> 7c, the two fixed points meet at A 0 . Finally, when 7 < 7 C , the zeroes move off 
the real axis, and fixed points no longer occur. Note that near the IR fixed point, 7 ~ e, 
whereas 7 C ~ y/e/N. Thus, for sufficiently large N, the critical value of the anomalous 
dimension needed to have finite BCS fixed points can be made arbitrarily small. Therefore, 
in this theory, N acts as a tuning parameter in the space of theories for 7 C . 

Next, consider what happens as 7 —» y c from above. In this limit, the UV and IR fixed 
points approach one another, and when this ratio becomes unity, the fixed points annihilate. 
Once this happens, for 7 < y c the system develops a superconducting instability, and the 
metallic phase is lost. In this case, the fermion anomalous dimension is not strong enough 
to avoid a superconducting instability. The inverse correlation length associated with the 
BCS coupling can be estimated as follows (see also Ref. 34 where such behavior is studied in 
detail): 


4c - A exp 


r r x " d\ 


A™ 

~ exp 


7T 


2\/7c T 2 


A 


(IV. 3) 


The correlation length diverges exponentially as 7 c /7 —> 1, signaling a continuous phase 
transition that separates the quantum critical and superconducting states. This behavior is 
similar to the way in which the correlation length of the 2d XY model diverges at the BKT 
transition. The analog of critical temperature is played here by 7 C . 


Notice that for 7 0, (IV.3) reproduces the boson-enhanced scale (IV.1). On the 

other hand, as 7 ~ 7 C , we find strong non-Fermi liquid corrections to the superconducting 
order parameter, which eventually destroy it via a continuous phase transition. Passing 
this phase transition obtains a quantum critical point characterized by non-Fermi liquid 
exponents for the quasiparticle dimension, Yukawa coupling, and BCS interaction. This 
QCP is characterized by a finite BCS coupling, with no superconductivity, and with a 
power-law behavior for the superconducting correlation function. 
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B. Quantum criticality and superconductivity 


Before, we presented a qualitative discussion of Eq. (III. 11). Our task now will be to 


study in more detail the phase structure and low energy dynamics as a function of N and e 
To begin with, let us start from the non-Fermi liquid fixed point at a = 2y = e/3. Then, 
the discriminant of the BCS j3 function vanishes at 


eN = 12. 


For N > 12/e, the BCS coupling flows to the stable IR fixed point 

A_ = yeiV(l - v/1 - 12/(dV)). 


(IV.4) 


(IV.5) 


Note that A_ ~ 0(1) over all the critical range. An important point here is that the 
attractive fixed point has a finite domain of attraction: for sufficiently large initial values of 
A, the suppression from the anomalous dimension term does not set in fast enough and A 
will diverge in the IR, signaling a superconducting instability. This is due to the existence 
of the unstable fixed point at A + . The size of the domain of attraction at the NFL fixed 
point is of order |A+ — A_|. This is much bigger than 1 for N 12/e, and shrinks to zero as 
the continuous transition is approached at IV = 12/e. These results are explained in more 
detail in the Appendix. 

It is possible to solve exactly the coupled system of equations 


do. e n 

= ",« + « 

dfi 3 


A 2 


d\ , 

^ = “ 2,r “ + “ A “ 2tt 2 N 


(IV.6) 


in terms of hypergeometric functions. The boundary conditions for the RG are imposed at 
some high scale n — M which should be below Mo given our approximations on the Zb — 3 
boson scaling. In accordance with our previous analysis, this obtains a fixed point a = e/3 
and A = A_ for N > 12/e, a superconducting phase for N < 12/e, and a continuous phase 
transition at N = 12/e. 

Near the continuous phase transition at IVe = 12, the scale of the superconducting insta¬ 
bility (namely the scale at which A diverges) is found to be 


hsc — exp 


7r /31V 
2 


y/1 - elV/12 


A. 


(IV.7) 
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For eN -C 12 this shows the enhancement due to boson exchange; however, as eN —>■ 12 
non-Fermi liquid effects dominate over this enhancement and destroy the superconducting 
parameter by a characteristic BKT scaling. This is one of our main results regarding the 
competition between superconducting and non-Fermi liquid effects due to a critical boson. 

An illustrative way of presenting these RG results is in terms of streamlines for (/ 3\, f3 a ) 
as a function of (A, a). The left panel in Fig. [b] shows the case N > 12/e, and the red 
points correspond to the stable, unstable and Gaussian fixed points. The flows in the 
superconducting range N < 12/e are presented in the right panel, where we also show the 
Gaussian fixed point. 


SJ 


A 

FIG. 6. Streamlines for (/3\,f3 a ), with N > 12/e (left), and N < 12/e (right). The red dots are 
fixed points; we show the stable, unstable and Gaussian fixed points in the left panel, and the 
Gaussian fixed point in the right panel. We thank A. Maharaj for help in generating this plot. 




V. DISCUSSION 

In this paper, we have studied a class of quantum metals (such as the Ising nematic 
quantum phase transition), obtained by coupling a Fermi surface to a nearly critical bosonic 
order parameter which preserves translation invariance, and condenses at zero momentum. 
We analyzed the interplay between superconductivity and non-Fermi liquid effects in a 
theoretically controlled setup depending on N (the rank of an internal global symmetry) and 
e (the deviation from spatial dimension d = 3). We found a novel class of fixed points, stable 
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against the superconducting instability, where the BCS interactions flow to scale invariant 
values. These QCPs display non-Fermi liquid behavior in observables that are local on 
the Fermi surface (the anomalous dimension and Fermi velocity) but also in operators that 
combine fermions on antipodal points, such as the Cooper pair field or the BCS operator. 
We also showed that for sufficiently small N a superconducting instability sets in, via a 
continuous phase transition. We next consider the possible relevance of our findings to the 
experimental observations of superconducting domes near quantum critical points in a broad 
class of correlated electron materials. 




(a) (b) 


FIG. 7. Fate of the superconducting domes for different values of Ne: (a) Ne <C 1 (left), and (b) 
Ne > 12. Here the various labels denote the following: the ordered phase (I), Fermi liquid (II), 
and non-Fermi liquid (III). 

We consider the phase diagram as a function of temperature and the parameter that 
tunes the boson to criticality (e.g. doping, pressure), which we label x. For N <C 1/e, the 
massless boson produces a strong enhancement in the BCS interaction, but a negligible non- 
Fermi liquid anomalous dimension. In this case we find a superconducting dome covering 
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the critical point. As N is increased, non-Fermi liquid effects become stronger, with the 
result that the fermion anomalous dimension tends to make the 4-Fermi attraction irrelevant 
and decreases the scale of superconductivity. Finally, for N > 12/e, superconductivity is 
destroyed and the non-Fermi liquid fixed point emerges. The competition between non- 


Fermi liquid and superconducting fluctuations is summarized in the scaling (IV. 7) for the 
gap. These different regimes are presented schematically in Fig. [7} 

In this work, we have considered a large N theory where the fermion is in the funda¬ 
mental of SU(N) and the scalar is in the adjoint. The N x N bosonic degrees of freedom 
were crucial for obtaining a quantum critical point stable against the BCS interaction. An¬ 
other possibility, discussed before in e.gJ 16 ^ ^ 0|33|35 | 3ti l , is to introduce a large number Np 
of fermionic fields while the order parameter remains a singlet. In such theories, it seems 
unlikely to obtain superconducting domes that condense out of a non-Fermi liquid normal 
state. The primary reason for this is that the fermion anomalous dimension in such theories 
is proportional to a/Np. By contrast, in our formulation, the fermion anomalous dimension 
7 ~ a. Comparing scales, in large N F theories, 


^NFL ~ e 


A 


(V.l) 


whereas the scale at which the superconducting instability develops is the same as in our 
theory, 


_ 7T j Np 

/4c ~ e 2 V - A . 


(V.2) 


It therefore is virtually impossible in a perturbative framework for non-Fermi liquid behavior 
to occur at scales above the superconducting instabilities in the large N limit of this class 
of theories. A similar conclusion, though from a somewhat different approach, has been 
reported in 33 . 

In future work, we wish to study the phenomenological consequences of the IR stable 
fixed point with finite BCS coupling. In particular, it would be interesting to determine 
the effect of a scale-invariant BCS interaction on thermodynamics (e.g. heat capacity) and 
transport (e.g. magnetoresistance, resistivity) properties. We wish also to study the effect 
of a magnetic held in our scenario. For instance, in the regime of N, e where there is a 
superconducting dome enveloping a quantum critical point, it is natural to ask what the 
properties of the system are when a magnetic held is used to destroy superconductivity. It 
is conceivable that when the superconducting dome is destroyed by a magnetic held, there 
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still remain substantial superconducting fluctuations governed by the BCS coupling; we can 
then ask whether the Cooper pair fields retain power-law correlations. 

We also found regimes where even though Ne > 12, for sufficiently small initial a or for 
large enough A, the system still flows towards the superconducting phase; see Fig. [6j This 
gives to a superconducting instability deep inside the non-Fermi liquid state. It would be 
interesting if this is of relevance to the phase diagram of some of the high T c superconductors. 
Finally, the analysis of the superconducting phase using the gap equation, and its connection 
with the present RG approach, will be discussed in 30 . 
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Appendix A: Scaling analysis 

In this Appendix we construct the scaling theory (within the spherical RG) that agrees 
with the form of the SD quantum corrections. 

In the presence of a nontrivial dynamical exponent Zb, the RG approach of Ref. ESI needs 
to be modified. First, the scaling of the fermions is not changed by zy from the action near 
the Fermi surface, 

Sf = — J dpodp±d d ~ 1 hi/)* r (ipo — vp±)'ijj , (A.l) 

we read off the scaling dimensions 

bo] = bJ = 1 , [b(p)] = -3/2 . (A.2) 

On the other hand, the bosonic momenta that dominate quantum corrections appear as 
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differences of close-by fermionic momenta, as can be seen from the cubic interaction, 

SVuk = J dpodp ' 0 d d p d d p g <i>{;p : - p)^ {p')f(p) , (A.3) 

Given p = h(k F + p± ), let us decompose the other fermion momentum as 

p' = h'(k F + p'±) ~ h{k F + p' ± ) + k F dh . (A.4) 


The boson momentum decomposed with respect to the local Fermi surface direction n then 
satisfies 

q = hq ± + gj| , q ± = p ± - p' ± , gj| = k F Sh . (A.5) 


Therefore, [g 0 ] = [q±] = 1, and it remains to understand how to scale 5h. 


The scaling of gj| is determined by the Zb = 3 boson propagator (II.10). Since [q±] = 1, it 
is gy that is affected by the dynamical exponent, and hence is the component that dominates 
the momentum transfer. We conclude that for bosonic momenta, 

10-e 


[go] = [q±] = 1 , [q\\] = 1/3 , 


6 


(A.6) 


for d = 3 — e. The scaling of the bosonic momenta has become anisotropic due to the 
dynamical exponent. 

From a purely bosonic point of view it seems somewhat artificial to select a direction 
n and scale the two components q± and q\\ differently, as this breaks the isotropy of the 
boson dispersion relation. Our statement, however, is that this is the scaling that will 
dominate inside correlation functions, where the boson momentum behaves as a difference 
of two fermionic momenta. In many condensed matter systems, the bosons indeed represent 
fermionic collective modes, hence it is very natural to identify their momenta with the 
difference between those of fermions. The scaling of q» is then equivalent to scaling differences 
in fermionic angles k F Sh, a process that determines the size of the Fermi surface patch that 
couples more relevantly to a given fermion ip(p,n). We note the related RG analysis of the 
fermion-boson system in the patch picture 1 -^, though we stress that the spherical RG being 
used here is not the same as the patch scaling of 1735 36 . 

With these scalings, the classical dimension of g calculated from 

Sb—f dp 0 dp' 0 dp±dp' ± d 2 ~ e [n + h')d 2 ~ e (h — h') g <f{p — p)f>\p)i>{p) • ( A -7) 


becomes 



(A.8) 
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Appendix B: RG solution for the BCS coupling 


We found that the stable fixed point A_ has a finite domain of attraction due to the 
existence of the unstable fixed point at A+. One interesting consequence of this is that we 
could have an RG trajectory that ends in a superconducting instability even if N > 12/e. 
This may then realize a superconducting dome condensing out of a non-Fermi liquid and 
could be of relevance for certain strongly correlated materials. For this reason, in this 
Appendix we discuss in more detail how this occurs in a simple case. 

By fixing a to its critical value a = e/3, we can solve the RG analytically across the 
transition, starting from arbitrary UV boundary condition A(M) = Ao- Denoting // = Me _t , 
the solution takes the form: 


A (t) = n 2 aN 
= n 2 aN 

= n 2 aN 


1 a 2 

a 

1 T 

~ N 

fo? 

a 

' T 

~ N 

/ 

a 

a 2 

I 

N ~ 



1 

2 

-] 

2 



aN > 4, A_ < Ao < A+ 
aN > 4, otherwise 

\ aN < 4. (B.l) 


From these solutions we can see that for aN > 4, the only possible pole occurs when the 
argument of the coth function becomes zero. Since the log term in the argument is positive 
for Ao < A_, for t > 0 a pole can only exist for Ao > A+. The corresponding superconducting 
scale is at: 

A = M (Z/A/) (B.2) 

Viewing Ao as a tuning parameter, there is a phase transition at A c = A+, whose order 

is given by , 1 . The order diverges as we tune aN towards the critical value 4. For 

V a2 ~w 

aN < 4, the superconducting scale is given by: 


where 


fi sc = Mex P 


2tt 


— - a 2 

N u 


/(Ao), 


/(A 0 ) = exp 


/ n 2 ^AaN — a 2 V 2 A 
\ n 2 aN — Aq / 


i / /a 

V N 4 


(B.3) 


(B.4) 
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We see that once we cross below aN = 4, the phase transition in Ao reduces to a step 
discontinuity across A disc — vr 2 aiV: there is no choice of Ao that can kill the superconducting 
instabilities. 
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